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Ch ■ Abstract 



We compute a number of distance-dependent universal scaling functions char- 
O ■ acterizing the distance statistics of large maps of genus one. In particular, 

we obtain explicitly the probability distribution for the length of the shortest 
non-contractible loop passing via a random point in the map, and that for the 
distance between two random points. Our results are derived in the context of 
bipartite toroidal quadrangulations, using their coding by well-labeled 1-trees, 
which are maps of genus one with a single face and appropriate integer vertex 
labels. Within this framework, the distributions above are simply obtained as 
scaling limits of appropriate generating functions for well-labeled 1-trees, all 
expressible in terms of a small number of basic scaling functions for well-labeled 



^ ■ plane trees. 

(N 
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O ' 1. Introduction 

2 ■ The understanding of random maps is a fundamental issue in combinatorics and 

O . many map enumeration results were obtained over the years, using for instance recursive 
decomposition |]T|, matrix integrals ^ or bijective methods [^0. Recall that a map 
of genus /i is a proper embedding of a graph in Sh, the compact oriented surface of genus 
^ ' h without boundary. By random maps, we mean generically some statistical ensemble 
■ of maps distributed according to some particular given law. Of particular interest in the 
scaling limit of large random maps, which converge toward nice universal probabilistic 
objects whose metric properties are only partially understood. This limit is especially 
relevant to physics, where maps are used as discrete models for fluctuating surfaces, e.g. 
in the field of biological membranes or that of string theory. In physics, maps are often 
equipped with statistical models, such as spins or particles, presenting a large variety of 
critical phenomena. This gives rise to many possible sensible scaling limits of continuous 
surfaces, each defining a particular universality class of maps. Their classification and 
characterization are the aim of the so-called two-dimensional quantum gravity [^,0, 
which also has deep connections with SLE processes P,p!0|. 



The simplest universality class is that of the so-called pure gravity, which describes 
the scaling limit of large maps with prescribed face degrees, such as triangulations 
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(maps with faces of degree 3 only) or quadrangulations (maps with faces of degree 4 
only), possibly equipped with non-critical statistical models. In this universality class, 
a particular attention was payed to the limit of large planar maps, i.e. maps with 
the topology of the two-dimensional sphere, converging to the so-called Brownian map 

Even if its topology remains that of a sphere [OJn], the Brownian map presents 
nevertheless intriguing metric properties such as a remarkable phenomenon of confluence 



of its geodesies [15|, which reveals an underlying tree-like structure. More quantitatively. 



the geometry of the Brownian map was characterized by a number of distributions 



measuring its distance statistics, such as its two-point function |16-19], which is the law 



for the distance between two random points in the map and its three-point function , 
which is the joint law for the three mutual distances between three points. Other reflned 
distributions, measuring e.g. the length of "separating loops" or that of the common 



part of confluent geodesies, were also obtained [^. Beside the spherical topology, the 



case of maps with a single boundary was also considered , with an explicit derivation 
of the law for the distance of a random point to the boundary. 

The most advanced results on the distance statistics in maps were obtained in the 
context of planar quadrangulations, using a bijection by Schaeffer [|| which gives a cod- 
ing for these maps by so-called well-labeled trees, i.e. plane trees whose vertices carry 
integer labels with particular constraints. This approach turned out to be the most 
adapted to address distance-related question since the labels precisely encode some of 
the distances in the original map. An extension of the Schaeffer bijection was also 
found by Marcus and Schaeffer |2^ - |2^ , which establishes a bijection between bipartite 
quadrangulations of arbitrary genus h and so-called well-labeled h-trees, which are maps 
of genus h with one face only, and whose vertices carry integer labels with the same 
constraints as in the planar case. In view of this result, a very natural question is there- 
fore that of the distance statistics in maps of arbitrary genus. The Marcus- Schaeffer 
bijection was used very recently [|2^ to discuss the scaling limit of large bipartite quad- 
rangulations of flxed genus, showing in particular the convergence to a limiting object 
with Hausdorff dimension 4, yet to be characterized. 

The purpose of this paper is precisely to give a quantitative characterization of 
the distance statistics in large toroidal maps, i.e. maps of genus 1. Here it is worth 
mentioning that toroidal cells or vesicles do exist experimentally |j2^ and some of our 
results might be relevant to their description. We will in particular derive an explicit 
expression for the two-point function of large toroidal maps, which differs from that of 
the planar case. In the case of genus 1, we may also deflne a non-trivial "one-point 
function" by measuring the length of the shortest non-contractible loop passing via a 
given point. We give here the corresponding limiting probability distribution for large 
maps. These results are obtained again in the simplest context of genus 1 bipartite 
quadrangulations, using the Marcus-Schaeffer bijection above with well-labeled 1-trees. 

The paper is organized as follows: in Section 2, we present our basic tools which 
are a number generating functions for well-labeled (plane) trees with marked vertices, 
with a particular emphasis on their scaling limit. We then recall in Section 3 the 
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Marcus-Schaeffer bijection between pointed (i.e. with a marked vertex called the origin) 
bipartite quadrangulations of genus h and well-labeled /i-trees. As an exercise, we 
show how to recover the number of genus 1 bipartite quadrangulations from our basic 
generating functions of Section 2. We then compute in Section 4 two different one-point 
functions for large pointed toroidal maps: the probability distribution for the length of 
the shortest non-contractible loop passing via the origin of the map and that for the 
length of the " second-shortest" loop, which is the shortest among non-contractible loops 
not homotopic to the shortest non-contractible one (or its powers). We finally derive 
in Section 5 the expression for the two-point function of large toroidal maps, which is 
the law for the distance from the origin of a random point in the map. We gather a few 
concluding remarks in Section 6. 



2. Basic tools: generating functions for well-labeled trees 

In this section, we recall a number of known expressions for generating functions 
of well-labeled trees and derive a few new ones. As will be apparent later, the explicit 
formulas displayed here will serve as basic tools to obtain, in the next sections, various 
distance related generating functions in ensembles of large toroidal maps. 

By well-labeled tree, we mean a plane tree whose vertices carry integer labels such 

that: 

(i) the labels on two vertices adjacent in the tree differ by at most 1; 

(ii) the minimum label is 1. 

It is useful to also introduce almost well-labeled trees where (ii) is replaced by the 
weaker condition: 

(ii)' all labels are larger than or equal to 1. 

The first generating function of interest is that, Ri = Ri{g), of planted almost well- 
labeled trees, with a weight g per edge, and with root label I {(->!). In other words, 
the coefficient of g'^ in the expansion of i?^ as a power series in g gives the number of 
planted almost well-labeled trees with n edges and root label I. This generating function 



was computed in |T^, with the result: 

(l-X^)(l-X^+3) 



Rp — R 



(l-X^+l)(l-X^+2) 



with R = and a; -I h 1 ^ 



(2.1) 



Qg X gR"^ 

This particular form is obtained by solving of the recursion relation 

with initial condition Rq = 0, expressing that an almost well-labeled tree with root label 
£ is coded by the sequence of its descending subtrees, themselves almost well-labeled 
trees with root label £ or 



3 



The contribution to Ri of large trees is encoded in the behavior of this generating 
function in the vicinity of the critical point g = 1/12. Setting 



^7 = ^(1-6^) (2.3) 
with e — 7- 0, a non-trivial scaling limit is reached for Ri by letting i become large as 

(2.4) 

with L finite. In this limit, we have the expansion 

Re ^2 (l-eJ^{L) + 0{e^/^] 

3 (2.5) 



with J'(L) = 1 + - 

sinh^ U^L 

The scaling function J-" satisfies the non-linear differential equation 

jr" = 3(jr2 _ 1) (2.6) 

which is the continuous counterpart of eq. (|2.2| ), obtained by expanding this discrete 
equation at order e^. 

Another basic generating function is that, Xi {£ > 1), of almost well-labeled trees 
with two distinct (and distinguished) marked vertices carrying the same label £ and such 
that the labels on the (unique) shortest path in the tree joining these two vertices are 
all larger than or equal to i. Adding for convenience a trivial contribution 1 to Xi, we 
have the recursion relation 

Xi = l + gRjXiil + gRj+^Xi+i) (2.7) 

obtained by inspecting the first occurrence of a label i on the shortest path joining the 
two marked vertices. The solution of (|2.7| ) reads 



(l-x3)(l-x^+i)^(l-x^^+3) 

^~ (l-x)(l-X^+3)2(l-x2^+l) ^ 



with X as above. In the scaling limit, we have the expansion 

Xt = 3(l--/~eC(L)+0(e)) 



with C(L) = 76 ^±^^2!^ 
sinh(v6L) 



while eq. (|2.7|), expanded at leading order in e, translates into 

C' = C2-6J^. (2.10) 
Prom the explicit forms of T and C above, we observe that me may write the identity 

C = -^ = -(log(-r))' (2.11) 



since \T'{V)\ = —J-''{L) for real L. In this form, the relation (|2.10|) is a direct conse- 
quence of the relation ( [2.61 ). 

A third useful generating function is that, X^i,^2 (^i > ^2 > 1), of almost well- 
labeled trees with two distinct (and distinguished) marked vertices carrying respective 
labels ii and and such that labels on the (unique) shortest path in the tree joining 
them are all strictly larger than £2 (except of course at the endpoint of the path with 
label £2)- By cutting this shortest path at the first occurrence of a label £1 — 1, £1 — 2, ■ ■ •, 
we may write 



ii-i 

£=£2 

(1- a;^^)(l- a:^^+i)(l- x^^+^)(l- 3;^^+3)(l- x^^i+3) 
~ ^ ' (1- a;^i)(l- a;^i+i)(l- x^i+2)(l- a;^i+3)(l- a;2^2+3) 



:2.12) 



This last expression is extended to the case £1 = £2 by adopting the convention that 
-^f-2,£2 — 1- scaling limit, setting £1 = Li/y^ and £2 = L2/y/e, we have the 

expansion 

Xi„i,=C{Li,L2) + 0{^e) 



cosh (^'^L'^ sinh^ (\/f -^1 

Note that this expression can be obtained alternatively by looking at the continuous 
counterpart of the product formula in ( p.l2|) . Indeed, using gRiRi^iXi^i ~ (1/12) x 
2^ X 3 (1 — i/eC(L)), we may write directly 

C(Li,L2) =exp l\{L)dL^ = exp (- [ - log(-^')] = (2-14) 

which matches (|2.13| ). 

The final important quantity is what we shall call the "propagator", which is the 
generating function Ki-^j^ (^1, ^2 > 1) of almost well-labeled trees with two distinct (and 
distinguished) marked vertices with respective labels £1 and £2, and with no additional 



condition on the labels in-between (apart of course from the general conditions (i) and 
(ii)' of almost well-labeled trees). This generating function reads 

min(£i,£2) 

Kii,t2 = -^£1,(2 + ^ei,e^£2,£^£ (2-15) 

£=1 

as obtained by summing over the minimum label £ on the path joining the two marked 
vertices. Note that we insist here on the two marked vertices being distinct, with a path 
of non-zero length in-between and that, if ii = £2, we have to subtract an undesired 
contribution 1 in the sum. The generating function Kg-^^g^ is solution of the equation 

^£1/2 = 9R£i{R£i+liS£i+l,£2 +-^^1+1,^2) + R£iiS£i,£2 +-^^1,^2) 

X (2-16) 

+ -R^l-l (5^1 -1,^2 + -^^1-1,^2 



with Ko^e^ = 0, obtained by looking at the label £1, £1 ± 1 of the vertex adjacent to the 
first marked vertex on the path joining the two marked vertices. We have no compact 
formula for Kg-^^^i^ at the discrete level but, in the scaling limit, eq. (|2.15| ) translates 
into the expansion 

p(Li,L2) ^ 
J^£i,£2 = -7^ ^ ^(IJ 

■ <T T ^ (2-17) 



with p(Li,L2) = / 3C{LuL)C{L2,L) dL 
Jo 

with a simple factor 3 for the limit of Xg. Using ( |2.14D , we may now compute the 
integral and write a compact expression for p(Li,L2): 

p(Li, L2) = r (Li)r (L2) X min(Li, L2)) 

where nil) = ^SOL + VQ ^sinh (^2y6L) -16sinh (^Vgl) -32tanh ^^f^ 

'(2.18) 

Note that the somewhat involved function 'H(-L) is simply characterized by 

n'iL) = ^ (2.19) 

with ?^(0) = 0. It is a straightforward exercise to check that the continuous propagator 
p(Li, L2) above is solution of the equation 

^p(Li, L2) = 6^(Li)p(Li, L2) - 3 S{L, - L2) (2.20) 



which is the continuous counterpart of (|2.16| ). Up to a (somewhat arbitrary) global 
normalization of the coordinate L, this is precisely the equation satisfied by the two- 
point function of a polymer chain embedded in one dimension (with coordinate L) 
subject to a potential V{L) oc J-'{L) — 1 oc 1/ sinh^(A/3/2L). Heuristically, this polymer 
corresponds to the (unique) path between the two marked points in the tree and its 
configuration in one- dimensional space simply reproduces the sequence of labels along 
this path. From condition (ii)', the presence of the other attached subtrees exerts on the 
polymer an effective repulsion from the position L = encoded in the potential V{L). 

When £i oo, Ki^^i^ — )• unless we also send I2 00, keeping p = £2 — ^1 finite. 
Defining 

kp = lim Ki^^i^+p , (2.21) 
for an arbitrary integer p, this quantity satisfies the equation 

kp = g R^{kp-i + kp + kp+i + 5p,i + Sp^ + dp-i) (2.22) 
and we find explicitly 

= TT^W ^ - ^P^^ (2.23) 

(1 — x)(l — X^J 

with X as in ( |2.1| ). In the scaling limit, and for large values of Li and L2, we have the 
corresponding limiting behavior 

/r 

p{Li,L2) ~ ^exp(-V6|Li -L2I) (2.24) 

where we recognize the usual (translation invariant) two-point function for a polymer 
in one dimension without potential. 

The above quantities Xi, X^^/j K^-^^i^ may alternatively be viewed as generat- 
ing functions for chains of almost well-labeled trees, with endpoints carrying prescribed 
labels. In the following, we will consider various generating functions for more in- 
volved well-labeled structures coding for toroidal maps. At the discrete level, those are 
naturally formed of a number of chains as above, attached by their endpoints. The 
knowledge of the scaling forms C(L), C{Li,L2) and p{Li, L2) for the chain generating 
functions will be sufficient to derive explicit expressions for the scaling form of these 
more involved toroidal generating functions. This eventually will translate into explicit 
universal scaling functions characterizing the distance statistics of large toroidal maps. 



3. The Marcus-SchaefFer bijection for bipartite quadr angulations of genus 1 

Our approach on the distance statistics of large toroidal maps relies on the bijection 
by Marcus and Schaeffer between, one the one hand, pointed bipartite quadrangulations 
of fixed genus h and on the other hand, well-labeled h-trees. By pointed bipartite quad- 
rangulation, we mean a map whose faces all have degree 4 (quadrangulation) , whose 
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vertices can be colored in black and white in such a way that adjacent vertices have a 
different color (bipartite) and with a marked vertex called the origin (pointed). It was 



shown by Marcus and Schaeffer ||2^ that any such map of genus h is bijectively coded 
by a well-labeled h-tree, i.e. a map of genus h with exactly one face [h-tree) and whose 
vertices carry integer labels subject to the two conditions (i) and (ii) of section 2 (well- 
labeled). As before, it is convenient to also introduce the notion of almost well-labeled 
h-tree where condition (ii) is replaced by the weaker condition (ii)'. 

Note that a quadr angulation with n faces and genus h has 2n edges and n + 2 — 2h 
vertices. The n faces of the quadr angulation are in one-to-one correspondence with the 
edges of the well-labeled h-tree, which therefore has n edges and n+l—2h vertices. These 
n-\-l—2h vertices are in one-to-one correspondence with the {n-\-2—2h) — l vertices of the 
quadrangulation other than its origin, and the label of a vertex in the h-tree is nothing 
but the distance from the associated vertex to the origin in the quadrangulation. Finally, 
the 2n corners of the h-tree (i.e. the angular sectors between consecutive edges around 
a vertex) are in one-to-one correspondence with the 2n edges of the quadrangulation. 
More precisely, the corners with label i (i.e. around a vertex with label £) in the h-tree are 
in one-to-one correspondence with the edges of type (^ — 1) — > £ in the quadrangulation, 
i.e. the edges connecting a vertex at distance {£ — 1) from the origin to a vertex at 
distance i. 

It is useful to recall how to recover the quadrangulation from its well-labeled h-tree 
coding. For the vertices of the quadrangulation, we take all the vertices of the h-tree 
plus an extra vertex added inside its unique face. This vertex will be the origin of the 
(pointed) quadrangulation. The edges of the quadrangulation are obtained by linking 
each corner with label i in the h-tree to its successor, which is the added vertex for £ = 1 
and, for i > 1, the first corner with label i — 1 encountered, say counterclockwise around 
the unique face of the h-tree. These links can be drawn without mutual crossings inside 
the face of the h-tree and, by construction, they do not intersect the original edges of 
this h-tree. Note that the sequence of links between the successive successors of a given 
corner with label i provides a particular geodesic (i.e. shortest) path of length i in the 
quadrangulation, leading from the vertex underlying this corner to the origin vertex. 
After drawing the quadrangulation edges, we may erase all the original edges of the 
h-tree as well as all the labels. Note that the obtained quadrangulation is automatically 
bipartite. 

In the spherical case h = 0, the well-labeled 0-trees are nothing but the well-labeled 
trees of Section 2, so that our generating functions above translate into generating func- 
tions for pointed quadr angulations (note that we need not specify that the quadrangu- 
lation is bipartite since this is automatic for h = 0). For instance, choosing a corner 
with label £ in a well-labeled tree amounts to picking an edge at distance i from the 
origin (i.e. an edge of type (£ — 1) — > i) in the associated quadrangulation. The gener- 
ating function for pointed planar quadrangulations with a marked edge at distance i is 
therefore that of planted well-labeled trees with root label i. It is given by Rg — Ri-i 
since we have to subtract from Rg those configurations with minimum label larger than 
or equal to 2. These undesired configurations correspond exactly to almost well-labeled 
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^1 




t2 




(b) 



Fig. 1: The two possible backbones for 1-trees: (a) a generic backbone 
with two vertices of degree 3 carrying positive labels l\ and linked by 
three edges and (b) a degenerate backbone with one vertex of degree 4 
with positive label £, linked to itself by two edges. 

configurations having all their labels shifted by 1, hence they arc counted by Ri-i- We 
will use the same shift argument in the following to transform generating functions for 
almost well-labeled objects into generating functions for their well-labeled analogs. For 
instance, the generating function for pointed planar quadr angulations with two marked 
distinct (and distinguished) vertices at distance £i and £2 from the origin is given by 
K£-^/^ — Ki-^-i^e^-i in terms of the propagator. 

In this paper, we focus on the case h = 1 of toroidal topology. Well-labeled and 
almost well-labeled 1-trees may be classified according to their backbone obtained as 
follows: we first delete recursively each vertex with degree 1 and its incident edge until 
no degree 1 vertex is left in the 1-tree. The resulting object is an almost well-labeled 
1-tree whose vertices all have degree larger than or equal to 2 by construction. We call 
this 1-tree the skeleton of the original 1-tree. We may then erase each vertex of degree 2 
in the skeleton and concatenate the two incident edges, leading to a labeled 1-tree with 
vertices of degree larger than or equal to 3 and with arbitrary positive vertex labels. We 
call this 1-tree the backbone of the original 1-tree. From the Euler relation for maps, it 
is easily seen that only two types of backbones are possible: 

(a) generic backbones, made of two 3-valent vertices with labels ^1 > ^2 > 1, linked by 
three edges (see figure l-(a)); 

(b) degenerate backbones, made of a unique 4-valent vertex with label i >1 linked to 
itself by two edges (see figure l-(b)). 

To obtain all the almost well-labeled 1-trees leading to a generic backbone, we simply 
have to replace each of the three edges of this backbone by an arbitrary non-empty 
chain of almost well-labeled trees with endpoints labeled by £1 and £2, as counted by 
the propagator Ki^ . The corresponding generating function reads therefore 



Here a factor 1/3 is necessary to avoid over-counting as, in the sum, we recover 3 
times the same configuration upon permuting cyclically the three edges of the skeleton. 




if £1 > £2 



if £i = £2. 



(3.1) 
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Similarly, an extra 1/2 factor is necessary to balance the possibility of recovering the 
same configuration by exchanging the two vertices of degree 3 if they have the same 
label. This simple over-counting argument does not hold for configurations presenting 
a symmetry and these configurations are therefore enumerated in (|3.1| ) with an inverse 
symmetry factor, which is a customary statistics. Note that the symmetric configura- 
tions are expected to become negligible for maps of large size n so that the symmetry 
factors have no effect on the scaling limit that we will discuss. 

If we wish to keep only the well-labeled 1-trees, we have to remove those configu- 
rations for which the minimum label is larger than or equal to 2 and, using the same 
shift argument as above, we now get a generating function 

r (3.2) 

with the convention Ki-^^ = Kq^q = 0. 

Summing first over the smallest label £2, then over the largest one ii, we get the 
total contribution Wi of all 1-trees having a generic backbone 
00 I" ^ ^ ii — i 

^2=1 




£1-1 \ ^ £1-2 



Q^£i,ii + 3 E ) ~ ( Q^ii-l,£i-l + 3 E ^ii-l 

a £iM q / j £i 




(3.3) 



6 '^''^ 3 

£2<£i 



1-^ 3 x3(l + 2x + 2x2 - 2x' 



g U . 3^ p 2(1-X)4(1+X)2 ■ 

Similarly, the total contribution W2 of all well-labeled 1-trees having a degenerate back- 
bone is 



^2 = El(^'iA-^'i-iA-i) 



4 

£1=1 



hm -K? . (3-4) 



A-^£i,£i 



1 2_ x2(l + 2x)2 



4 " 4(l-x)2(l + x)2 ■ 
Summing ( p.3|) and (|3.4| ), we get the generating function for well-labeled 1-trees which, 
from the Marcus-Schaeffer bijection is also that, Q^^\g) of pointed bipartite quadran- 
gulations of genus 1 

Q. -^l+^2-4^^_^^4(l+^)2 • (3-5) 
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Note that the symmetry factor of a weU-labeled 1-tree is also that of the associated 
pointed quadrangulation so that symmetric pointed quadrangulations are counted in 
qI^'* with their usual inverse symmetry factor. It is more customary to consider rooted 
maps, i.e. maps with a marked oriented edge, as they do not involve symmetry factors. 
The generating function for rooted bipartite quadrangulations of genus 1 is simply 

Qii^ = 4Qi') (3.6) 

since there are exactly twice as many edges as vertices in a genus 1 quadrangulation, 
each coming with two orientations. All the enumeration formulas above are consistent 
with those found in ref. P^ . 

If we are interested only in large maps, we can use the continuous analogs of ( |3.3|) 
and (^]5) giving the leading singularity of Wi and W2 , namely 



1 1 



Wi^^ lim ^/ dL2P^{Li,L2) 
Li^cx) 3 Jq 

^ 3^/1'^^ (x) 1^1'- 

W2^- hm -p^ Li, Li) = — \ — \ = — 
e Li^oo 4 4e \ 4 y SZe 

which can also be directly read off (|3.3|) and ( |3.4| ) with x ~ exp(— v^6e). Note that the 
degenerate case is less singular than the generic case so that the dominant singularity 
of qI^"* or Q^_}y^ comes from Wi only. 

Alternatively, we deduce from ( |3.7| ) the large n behavior of the number of well- 
labeled 1-trees with n edges having a generic or degenerate backbone 

12'^ 3 12" 

V^^ll,"-^, W2\g^r^—^. (3.8) 

96 32v7m 



At large n, the dominant contribution to qI^'* or Q^_^ comes from configurations with 
a generic backbone only and we recover the known asymptotics ||2^ 



1 9*^ 

Qa^|."=4Qi^t~^ (3-9) 
for the number of rooted bipartite quadrangulations of genus 1 with n faces. 



4. One-point functions for large toroidal maps 

For maps with a toroidal topology, we may define several interesting distance- 
dependent one-point functions. Starting with a pointed map, we may consider as a "self- 
distance" of the origin vertex the length of any shortest non-contractible loop in the map 
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Fig. 2: A non-contractible loop (thick red) drawn on a toroidal quadran- 
gulation must intersect the skeleton of the associated 1-tree at one of its 
vertices. 



containing this origin. In the case of bipartite quadr angulations, we denote this length 
by 11 as it is necessarily even. We may now use the coding of previous section and note 
that any non-contractible loop in the quadrangulation must intersect the skeleton of the 
associated 1-tree. Now the edges of the quadrangulation lie strictly inside the unique 
face of the 1-tree and do not cross its edges so that the intersection with the skeleton 
necessarily takes place at one of the skeleton vertices, say with label I' (see fig. 2 for 
an illustration). Since I' is the distance in the quadrangulation from this vertex to 
the origin, the length of the loop is necessarily larger than or equal to 21' so that the 
half-length £ of a shortest non-contractible loop is larger than or equal to the minimum 
label on the skeleton of the 1-tree. Picking now a vertex with minimal label Imm on 
this skeleton, this vertex has a degree larger than or equal to 2 by construction and we 
may consider two distinct corners around it. The two sequences of links between the 
successive successors of these corners define two paths of lengths £min leading to the 
origin, which do not intersect the skeleton and whose concatenation therefore creates a 
non-contractible loop of length 2£jnin- We deduce that 



£ = £^i„(Sk.) (4.1) 



where £jnin(Sk.) is the minimum label on the skeleton of the 1-tree coding for the pointed 
quadrangulation at hand. 

In the scaling limit, as already noticed, it is sufficient to consider the contribution 
of configurations leading to a generic backbone and we denote by ii = Li/y/e and 
£2 = the labels of its two vertices of degree 3, with ii > £2, and by mi = Mi/^/e, 

1712 = and 7713 = M^/ y/e the respective minimal labels on the three chains linking 

these vertices in the skeleton, with mi > m2 > ms so that £min(Sk.) = ms. With these 
notations, the generating function for almost well-labeled 1-tree whose skeleton has a 
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minimal label £jnin(Sk.) = I = Ly/e behaves when e — t- as 



M3 



dMi 



M2 



00 /•OO 

dL2 / dLi 3C(Li,Mi)C(L2,Mi) x 3 C(Li, M2)C(L2, M2) 

Ml J L2 



dMo 



'r{M2 



dMi 



X{L) 



with X(L) 



96 



Mo 



tanh 




x3C(Li,M3)C(L2,M3 

/•OO noo 
^Mi 



(4.2) 

Note, as in ( ^.IT] ), the factors 3 for the limit of each of the discrete terms X^^, 

and ■ Note also the factor 2 accounting for the two possible order of appearance of 

the minima mi, 1712 and ms when turning around say the vertex with larger label £1. 

To return to well-labeled configurations, we note that, as an alternative to the shift 
argument used in Section 3, we may equivalent ly use the property that almost well- 
labeled 1-trees whose skeleton has a minimal label i are in one-to-one correspondence 
with well-labeled 1-trees whose skeleton has a minimal label smaller than or equal to 
i. This correspondence is obtained by now shifting all labels in the almost well-labeled 
1-tree by a non-negative integer so that its global minimum label (in the whole 1-tree) 
becomes 1. From the Marcus- Schaeffer bijection, we deduce that ( [4 .21 ) may alternatively 
be viewed as the e — )■ leading behavior of the generating function for pointed bipartite 
quadrangulations with a shortest non-contractible loop passing via the origin of (rescaled) 
length smaller than 2L. 



a(r) 



~ r 




Fig. 3: Plots of the cumulative probability distribution a{r) and the asso- 
ciated probability density (j'{r) for the rescaled half-length r of the short- 
est non-contractible loop passing via the origin of a large pointed toroidal 
quadr angulation. 

From our one-point toroidal scaling function X(L), we may extract a (cumulative) 
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probability distribution in the ensemble of bipartite quadrangulations of genus 1 with 
fixed size n (= number of faces), in the hmit of large n. At the discrete level, fixing 
the size amounts to picking the coefficient g"' in the generating function at hand, which 
can be done by a contour integral in g. In the limit of large n, this contour integral 



translates into an integral over a real variable ^ upon setting (see |jT^ for details) 



9=^(^ + -] ■ (4-3) 



12 V n 

A sensible scaling limit is now obtained by rescaling i as 

i = rn^/^ (4.4) 

so that we can use in practice our continuous expressions above with e = — i^ / i/n and 
L = a/— i^ r. Normalizing by Q^^^\g'^, we get the probability 

a(r) 



- [ d^J,^eM-e) (x(V^r) -X(V^ 

TT Jo -K ^ 



n 



di 2n sin(v^r) sinh^(v^r) - sinh(v^r) sin3(v^r) '^^'^^ 

T^^^^ (cosh(v^r) - cos(v^r))4 



that the shortest non-contractible loop passing via the origin has a length smaller than 
2rn^/^ in the ensemble of pointed bipartite quadrangulations of genus 1 with fixed size 
n, in the limit n — oo. The cumulative probability distribution a{r) and the associated 
probability density a"'(r) are plotted in fig. 3. For small r, we have the expansion 

9r6 1431rio ,4, , 

It is worth mentioning that, in the continuous limit, the shortest non-contractible loop 
is expected to be unique, i.e. two shortest loops at the discrete level remain at a distance 
negligible with respect to n^^^ at large n. Moreover, as for geodesic paths in spherical 
maps, we expect a confluence phenomenon with the shortest loop made of an open part 
of (rescaled) length 2r — 26 and a common part of length 6 traveled back and forth in 
the vicinity of the origin. The joint law for 6 and r could be obtained in principle by a 
slight refinement of the above analysis. 

In the same spirit, for a pointed map of genus 1, we may first consider a shortest 
non-contractible loop as above and now study the length of any shortest loop among 
those non-contractible loops not homotopic to the shortest loop or to any of its powers. 
A similar quantity was introduced in an analyzed using heuristic scaling arguments. 



At the discrete level, this notion depends on the particular choice of the first shortest 
loop but we expect that this dependence is wiped out in the continuous limit in which 
the first shortest loop is essentially unique. We may as before assume that the well- 
labeled 1-tree coding for the quadrangulation has a generic backbone, i.e. has a skeleton 
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Fig. 4: Given a first non-contractible loop (thick red) intersecting only 
one of the three chains forming the skeleton, any other non-contractible 
loop non-homotopic to this first loop or to its powers (thick blue) must 
intersect (at least) one of the two other chains in the skeleton (possibly at 
an endpoint of the chain). 



with two 3-valent vertices linked by three chains. We may moreover assume that the 
minimum label on this skeleton is reached for a bivalent vertex, i.e. strictly inside one of 
the three chains. Indeed, all the other situations (degenerate backbone or minimal label 
at one of the two 3-valent vertices of the skeleton) correspond to degenerate cases which 
give sub-dominant contributions in the continuous limit. We now take for the first loop 
the concatenation of the two sequences of links between the successive successors of the 
two corners at the above bivalent vertex with minimal label. A non-contractible loop 
not homotopic to this loop or to its powers must necessarily intersect the skeleton at 
one of the two other chains so that its length is larger than twice the label of any vertex 
of the skeleton minus the first chain (see fig. 4 for an illustration). Taking a vertex with 
minimal label i^i^ in this set, the two sequences of successors of two corners at this 
vertex, once concatenated, create a suitable loop of length 2i^in- 

To summarize, denoting as before by Li and L2 the (rescaled) labels of the 3-valent 
vertices of the skeleton, with Li > L2, and by Mi, M2 and M3 the respective minimal 
labels on the three chains linking them, with Mi > M2 > M3, the length we are now 
interested in is nothing but the second minimum M2 so that, integrating over the other 
variables, we are now led to compute 

-/ dMs dMi dL2 rfLi 3C(Li,Mi)C(L2,Mi) X 3C(Li,M2)C(L2,M2) 

^ Jo J M2 J Ml J L2 

X 3C{Li,Ms)C{L2,Ms) 

M2=L 
3 



pLj q /"OO q /"OO /"OO 



= ^ with J{L) = -]-H{L)r{2L) 
with n as in (pT8|) . 



(4.7) 
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Fig. 5: Plots of the cumulative probability distribution a2{r) and the 
associated probability density (J2{r) for the rescaled half-length r of the 
"second shortest" non-contractible loop passing via the origin of a large 
pointed toroidal quadrangulation. 

Repeating the above arguments, we get the probability 



(T2 r 



96 

n 



Jo ^ 



(4.8) 



that, in the set of all non-contractible loops passing via the origin, the "second shortest" 
loop, i.e. the shortest among loops not homotopic to any power of the true shortest one, 
has a length smaller than 2rn^/^ in the ensemble of pointed bipartite quadrangulations of 
genus 1 with fixed size n, in the limit n — )■ oo. The cumulative probability distribution 
(72 (r) and the associated probability density o'2{r) are plotted in fig. 5. For small r, we 
have the expansion 

, . 11043 , 
+ C(r^^) . (4.9) 



(72 r 



It is interesting to comment on the small r behavior of both a{r) and <J2{f) which 
informs us on the regime of large but finite values of namely: 1 <ti i <ti n^/*. For such 
a finite value of i, the average number of vertices in a random quadrangulation giving 
rise to a shortest non-contractible loop of length smaller than 2i scales as 



na{r) 



n X r 



6 



n 



n 



1/4 



n 



(4.10) 



where we dropped the multiplicative prefactor. In particular, this quantity tends to 
as n~^/^ for large n. This is to be contrasted with, for instance, the average number 
(oc n(£/n^/'*)^ = i^) of vertices at a finite distance less than i from the origin of a 
pointed planar quadrangulation [|19| which remains finite at large n. The explanation 



for this phenomenon is that a prerequisite for the shortest non-contractible loop passing 
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Fig. 6: A schematic picture of a toroidal quadrangulation with (a) a small 
non-contractible loop or (b) two non-homotopic small non-contractible 
loops. The first case occurs with a probability decaying as n~^/^ and 
the second with a probability decaying as n~'^/^. 

via the origin to be finite is that the smallest non-contractible loop in the map be 
itself finite, while ( |4.10| ) essentially counts the vertices lying at a finite distance from 



this non-contractible smallest loop. The n dependence in ( [4.10|) is compatible with a 
probability of having a finite smallest non-contractible loop decaying as n~^/^ in the 
set of bipartite quadrangulations of genus 1 and fixed large size n (see fig. 6-(a) for an 
illustration). Similarly, from the small r behavior of o"2(r), the number of vertices giving 
rise to a second-shortest loop of finite length smaller than 2i behaves as 

/ £ \10 £10 

na2(r)~nxri° = n(^^J = . (4.11) 

We now interpret the n dependence in ( ^1.11|) as measuring, in the set of bipartite 
quadrangulations with genus 1 and fixed size n, the probability that a quadrangulation 
has its two non-homotopic smallest non-contractible loops finite: this probability tends 
to as n~^l'^ (see fig. 6-(b) for an illustration). 

These scaling behaviors are not surprising since, in case (a) of fig. 6, the quadran- 
gulation may essentially be viewed as a planar quadrangulation with two marked points 
which are glued together to create a handle, while in case (b), the quadrangulation is 
essentially a planar quadrangulation with one marked point. Recall that the number of 
planar quadrangulations of size n with a marked vertex is [|1| 

— \jlL - (4 12) 

2nn + l v^nV2 v • ; 

while that with two marked vertices is (asymptotically) n times bigger. Dividing by the 
number of total number (5^)|gn/(4n) ~ 12"/(96n) of quadrangulations of genus 1 gives 
a ratio scaling precisely as n~^/^ for two marked points and n~^l'^ for one. 
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To summarize, we deduce a contrario that a generic bipartite quadrangulation of 
gemis 1 has both its smallest cycles of order n^/^ so that its topology remains that of a 
genus 1 surface in the scaling limit. 

5. Two-point function for large toroidal maps 




Fig. 7: The "diagram" enumerating generic configurations contributing to 
the dominant singularity of Z^. Each edge of the diagram must be replaced 
by a propagator K with the indices of its endpoints. We must then sum 
over £i, I2 and £3, keeping (. fixed. 

The most natural measure of the distance statistics is via the two-point function 
which measures the density of points at a fixed distance from a given origin in the map. 
More precisely, we may decide to enumerate pointed bipartite quadrangulations of genus 
1 with now a marked vertex at distance £ (i > 1) from the origin. From the Marcus- 
Schaeffer bijection, this amounts to enumerating well-labeled 1-trees with a marked 
vertex with label £. As before, it is simpler to consider the generating function Zi for 
almost well-labeled 1-trees with a marked vertex with label £. From the Marcus-Schaeffer 
bijection, Zg is now the generating function for pointed bipartite quadrangulations of 
genus 1 with a marked vertex at distance less than or equal to I. 

The dominant singularity of Zi comes from quadrangulations leading to well-labeled 
1-trees with a generic backbone and, moreover, with a marked vertex with label £ lying 
outside their skeleton. The shortest path in the 1-tree from the marked vertex to 
the skeleton defines a chain whose endpoint on the skeleton has label, say £3, and is 
generically different from the two vertices of degree 3 in the skeleton, with labels £\ and 
£2- The dominant singularity of Zi is therefore that of the "diagram" of fig. 7, whose 
value is obtained by replacing each edge by the corresponding propagator K and by 
summing over £1, £2 and £3, keeping £ fixed, namely 

E E E K,,,,K,,,,,K,,,,, {K,,,,,f . (5.1) 

^1>1^2>1^3>1 
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Setting £ = i^/i/e and turning to rescaled variables in the sums, we get a dominant 
singularity 



2 /»CX3 /»00 l*00 

Ze — dLi / dL2 / dLs p{L, Ls)p{Li, Ls)p{L2, Ls) {p{Li, L2)f 

^ Jo Jo Jo (t- o\ 



e4 



with a new scaling function J^^^\L) which fully characterizes the two-point distance 
statistics in large toroidal maps. The above integrals can be performed explicitly and 
involve 12 sectors for the determination of the minimum in the formula ( |2.1(j| ) for p (2 
relative positions for L and L3 x 6 relative positions for Li, L2 and L3), reducible by 
symmetry Li -H- L2 to 6 sectors. After some cumbersome but straightforward calcula- 
tions, we find the relatively simple expression 

7^^\L) = Ao{L) + LAiiL) + L''A2{L) 

1 238 + 151 cosh(y6L) + cosh(2y6L) 



with AoiL) 
Ai{L) 
A2iL) 



sinhM^/fL 



5 100 sinh(\/6L) + 31 sinh(2y6L) (5.3) 
2048^6 sinh^ ^^l] 

75 3 + 2 cos 



1024 sinh^f^L 



This explicit formula constitutes the main result of this paper, with J-'^^^ playing for the 
toroidal topology the same role as J-" for the spherical case. It is worth noting that, from 
their explicit expressions, the three functions Aq, Ai and A2, satisfy the remarkable 
relation 

A'o - 4A[ + 2OA2 = (5.4) 
so that we may write in all generality the parametrization 

^0 = 4q!o , Ai = aQ + 5ai , A2 = ct'i . (5-5) 

In other words, the scaling function J-'^^^ may be written as 

J^^^\L)=4M{L) + LM'{L) 

(5.6) 

with M{L) = ao{L) + Lai{L) 
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in terms of a simpler function Ai. We have no explanation for this particular form nor 
for the meaning of M.{L) itself. From ( |5.3|) , we have explicitly 



ao{L) 



ai{L) 



1 238 + 151 cosh( + cosh(2v^L) 

(5.7) 



3072 sinh^r./lL 



25 /S^osh JfL 



V 2 ^.^1^5 //3^ 



At small L, although -^i and A2 diverge respectively as and a 

number of cancellations ensure that J^^^\L) vanishes as 

7-4 

jr(i)(n . (5.8) 

^ ^ 896 ^ ' 

We deduce from this formula that for a large, but finite i.e. in the so-called local limit, 
we have 

t , , „ , 12" 



. „ and therefore ZAgr, ~ £ (5.9) 

896e4 896e2 "^'^ 896 ^ ' 

at large n. Normalizing by the number Qi^^lg" ~ 12"^/96 of pointed quadrangulations, 

we find that the average number (V^)!^'* of vertices at a distance less than or equal to 
I of the origin in pointed bipartite quadrangulations of genus 1 behaves, for large but 
finite £, as 

This is exactly the result found for planar quadrangulations, i.e. in the spherical case 
/i = 0. As might be expected, the local limit is totally insensitive to the genus. The 
presence of handles may generically be felt only by traveling along distances of order 
n^/^ in the map. 

Alternatively, we may consider the scaling limit of pointed bipartite quadrangula- 
tions of genus 1 with fixed, large size n, but now with £ oc 71^/^. The probability <^^^\r) 
that a vertex picked uniformly at random in the quadrangulation be at a distance less 
than I = rn}l^ from the origin is given by 

^^jlI'>i r^i_e.p(_a(^(..(V3]?.)-^^^) (5.11) 

28 7r ^ '\ vv ; ggg j \ > 

where we singled out the contribution from the first term in the small r expansion of 
T^^\r), corresponding precisely to the local limit, as it is proportional to the improper 
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Fig. 8: Plots of the cumulative probability distribution ^^^\r) and the 
associated probability density $(^)'(r) for the rescaled distance r from a 
uniformly chosen random vertex to the origin of a large pointed toroidal 
quadrangulation. 




... ^yr) 




6 

r 



Fig. 9: Comparison of the probability densities $(^)'(r) and $('^)'(r) for 
the rescaled distance r from a uniformly chosen random vertex to the 
origin in large pointed toroidal and planar quadrangulations respectively. 
We plotted on the right the difference - 

integral (d^/^) exp(— ^^). Its value —in of this integral is dictated by the change 
of variable from g to ^. The integrals in the second and third lines of are now 

convergent integrals. Expanding further at small r gives 



^ ' 28 



15r 



10 



1242135r^'^ 



14560F 506970464v/^ 



+ 



(5.12) 



with, as discussed above, the same leading term oc r"* as for planar quadrangulations, but 
now a first negative correction of order r^^ instead of for the planar two-point function. 
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The two-point function therefore increases faster at smaU r for toroidal maps than 
for spherical ones. The probability distribution <^^^\r) and the associated probability 
density $(^)'(r) are plotted in fig. 8. A comparison with the corresponding genus 
probability density $'^'^)'(r), as computed in ref. [ffOfi is displayed in fig. 9. 



6. Conclusion 

In this paper, we have derived explicit expressions for a number of probability 
distributions characterizing the distance statistics of large toroidal maps. These distri- 
butions, obtained here in the context of bipartite quadrangulations, are expected to be 
universal (up to a non-universal global rescaling of r) and describe the distance statis- 
tics in more general ensembles of large toroidal maps in the universality class of pure 
gravity, such has maps with prescribed face degrees, possibly equipped with non-critical 
statistical models. 

Our main result is the explicit form (|5.3| ) for the two-point scaling function J^'^^\ 
which is the genus 1 counterpart of spherical two-point scaling function J-" of eq. (|2.5|) . 
In the same way as J-" satisfies the non-linear differential equation ( pj.6| ) , we may wonder 
if J^^^^ itself obeys some simple differential equation, possibly involving J-" as a source. 
Our method, which consisted in a direct computation of J^^^\ did not allow us to find 
such an equation. 

Among possible extensions of our work, let us mention the computation of the 
toroidal three-point function or more simply, that of a more involved two-point function 
now measuring the "second-shortest length between two points", i.e. the length of any 
shortest path among those paths not homotopic to the true geodesic. There seems to 
be no fundamental obstacle to the derivation of these laws but the calculations may 
rapidly become heavy. 

Another natural extension if of course that of maps with genus h > 1. Here a 
more fundamental obstacle occurs since, when h becomes large, we have to deal with a 
large number of diagrams, each involving a large number of propagators, which makes 
in practice our method unadapted. Alternatively, one may hope for the existence of a 
hierarchy of equations satisfied by the higher genus two-point functions, whose discovery 
would be a promising step in the quest for these universal scaling functions. 
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